function [IA_Network, consistent] = IA_Network_Add_Relation(IA_Network, Interval_Name_i, Interval_Name_j, IA_Relation)
%IA_Network_Add_Relation To Add R(i, j)
%   Detailed explanation goes here

[~, i] = IA_Network_Find_Interval(IA_Network, Interval_Name_i);
[~, j] = IA_Network_Find_Interval(IA_Network, Interval_Name_j);
%[IA_Network, consistent] = IA_Network_Path_Consistency(IA_Network, i, j, IA_Relation);
[IA_Network, consistent] = IA_Network_Path_Consistency_Hogge(IA_Network, i, j, IA_Relation);

% [IA_Network_hogge, consistent_hogge] = IA_Network_Path_Consistency_Hogge(IA_Network, i, j, IA_Relation);
% [IA_Network_matrix, consistent_matrix] = IA_Network_Path_Consistency_Matrix(IA_Network, i, j, IA_Relation);
% consistent = consistent_hogge & consistent_matrix;
% for n = 1:length(IA_Network_hogge.nodes)
%   for m = 1:length(IA_Network_hogge.nodes)
%     if ~IA_Relation_Is_Equal(IA_Network_hogge.relations(n, m), IA_Network_matrix.relations(n, m))
%       consistent = 0;
%       break;
%     end
%   end
% end
% IA_Network = IA_Network_matrix;
end